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I. SlMiARY OF RESULTS ACCOMPLISHED DURING PERIOD - MAY 16, 1975 - 
NOVEMBER 15, 1975 

A. Analytical Results 

The basic two year study plan is shown in TABLE I. The 
items indicated by asterisks were proposed for study in liie second 
(present) year. Ihe items indicated by a check marie have been 
ccmpleted during the six nonth reporting period. IVoo basic types 
of movable appendages have been considered throu^out the study 
1) a hinged type of fixed length vdiose orientation relative to the 
main part can change and 2) a telescoping type of varying length 
vhich is further subdivided into (a) the case vhere a large end 
ness is mounted at the end of an assuned nassless boom and (b) 
where the appendage consists of a uniformly distributed homo- 
geneous mass throu^out its length (Fig. 1). 

The equations of motion for the hinged system as previously 
developed have been linearized about the nominal equilibrium posi- 
tion vhere the boons are orthogonal to the nominal spin axis for 
the case of two dimensional motion (as simulated by Lang and 
Honeycutt^*) and tliree dimensional motion. Analytic stability 
erdteria are developed from the necessary condition on the sign of 
all the coefficients in the system characteristic equation. 


Please refer to references listed after Table I. 
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It is seen that, for the general case, for stability scm form 
of active danping must be present and that for limiting cases 
where the spin axis is an axis of synmetry, certain inequalities 
relating the hinge point offset coordinates to the moment of 
inertia ratio must be satisfied. For this same case it is seen 
that the system characteristic equation separates into two factors - 
a second order factor describing a mode where both hinged members 
move in phase as a unit, and the second factor represented by a 
more conrolex fourth order polynomial. 

lb date viscous danping at both hinge points has been included 
in the nathenatical model . Rate damping about the transverse axes 
will be included since the numerical results indicate that the 
coupling with only hinge dauping present does not result in favor- 
able nutation decay time constants. 

Equations of motion have been developed for the telescoping 
system where: (1) a single boom is linearly offset from the nominal 

spin axis and (2) two booms are linearly offset fixim the hub princi- 
pal axes - one boom from the spin axis and the second frcan one of 
the transverse axes. For the latter case, the two booms move at 
right angles to one anolher. (Figs. 2a, 2b) These equations of 
motion have been linearized about two desired equilibrium states: 

(a) a final spin rate about the ncminal spin axis only and (b) 
a state of zero inertial angular velocity. It is concluded that 
three-axis control with a single offset boom can not be achieved 
and that boom offset is required for any control of the linear 
system. 
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For a special case of a single boom offset from tiie hub axis 
of mass syirmetry, it is seen that the optiial control theory 
f«* the linear system using a quadratic performance index with 
unspecified final time can be ejomined analytically. For such a 
system liie jratrix Riccati equation can be solved in closed form 
and the control law gains obtained without the use of a computer. 
For this application (involving two axis control) the integrand 
in the cost functional is assumeci to be a weighted function of 
the transverse angular velocity plus the control. 

The problem of optimal control with a minimum time criterion 
has been examined analytically for the special case of a single 
offset boom vhere it is assumed that the initial conditions are 
such that the system can be driven to the equilibrium (rest) 
state with only a single switching naneuver in the bang-bang 
optimal sequence. For this system it is possible to obtain an 
analytical solution for the switching and final times in terms of 
the initial conditions and magnitude of the maximum value of the 
control force. 

B. Numerical Results 

The nonlinear equations of motion describing the deployment 
of the hinged system have been programmed for numerical integra- 
tion using the NOVA cerrputer at Reward University. A few typical 
results are illustrated in Figs. 3-6. Figs. 3(a) and 3(b) simulate 
the dynamic response of the system to an initial r^rturbation in 
Ihe hinge amgles of 0.1 radians with hinge dairping absent and then, 
present, respectively. 
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For the damped case a damping coefficient of 0.1 Ib/ft/sec is 
selected. The deployment of the system from the positiw vdiere 
the hinged members are initially parallel to the local vertical 
(«i=p) is simulated in Fig. 4(a) without damping, and Fig. 4(b) 

k 

with damping. The 'x' represents the maximum time simulated by 
Lang and Honeycutt. It is seen that without damping the hinged 
maiibers esdiibit a flapping- type motion as mcmentum is exchanged 
between the hinge and spin motions. With the damping coefficient 
selected the system can be fully deployed in about 10 seconds. 

An example of 1he three dimensional hinge-system dynamics is 
simulated in Fig. 5 (undamped) and 6 (with danping). Initial 
perturbations in both hinge angles and one of the transverse 
angular velocities are assumed. Allhou^ "the hinge danping is 
effective in reducing the amplitudes of the hinge motion, the 
time constants associated with the nutation angle decay (Fig. 6(b) ) 
are extremely long. It is clear that an additional form of nuta- 
tion damping must be added for effective removal of excess trans- 
verse rates. 

An illustrative numerical example, using the single offset 

boom as a wobble damper with the gains in the control based on 

the linear optimal control theory is shown in Fig. 7. The dotted 

curve shows the response of the nutation angle for such a system 

vdiere the physical parameters, mcments of inertias, masses 

and offset length are identical to the single control mass system 

6 7 

previously analyzed by Edwards and Kaplan. * 
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Ih their application a control law is selected such that the 
average rate of change of excess rotational kinetic energy is 
less than zero for eadi cycle. Using the same value of naximum 
amplitude (i.e. mtodmum boon length) it is seen that the nutation 
decay time constant can be ijjpDved by about an order of nagni- 
tude by using a control law based on linear optimal control theory. 


II. PRESD^TATION AND PUBLICATION OF RESULTS 

On November 6, 1975 an oral presentation was given at NASA- 
Langley to Mr. C. W. Martz, the Technical Monitor, Dr. W. Anderson, 

Head, Stability and Control Branch of the Flight Dynamics and Control 
Division, and to 6 other members of the Stability and Control Branch. 

Prior to this date Mr. Martz had been contacted several times by tele- 
phone to infom him of the progress . Copies of all visual aids used in 
the presentation (about 50 in total) were left at Langley. It is 
anticipated that another presentation will be given during the periixl 
December 1975 - January 1976 at which time a draft copy of a follow-on 
proposal could also be discxjssed. 

During the six month reporting period two papers based on last 
year’s work were accepted: 

(1) "Spacecraft Detuirbling Using ^fovable Telescoping 
Appendages," Peter M, Bainum andR. Sellappan, 

)(XV1 International Astronautical Congress, Lisbon, 

September 21-27, 1975, Paper No. 75-113. 

(2) "Dynamics of Spin Stabilized Spacecraft During Deploy- 
ment of Telescoping Appendages," by R. Sellappan and 
Peter M. Baimm, to be presented at the AIAA - 14th 
Aerospace Sciences Meeting, Washingjton, D.C., Jan. 26- 
28, 1976, Paper No. 76-185. 

In addition a paper based m a NASA-Grant during the 1973-74 period 
has been published: "Ihree Dimensional Motion and Stability of 

TWo Rotating Cable - Connected Bodies," by P.M. Bainum and K.S. Evans, 
Journal of Spacecraft and Rockets, Vol. 12, No. 4, April 1975, pp. 242-250. 
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III. PIANS FOR THE NEOT REPORTING PERIOD 

Additional nvinerical sijnulations of the hinged system 
during deployment and in response to snail perturbaticans will 
be made with both hinge danping and nutation damping present, 
and the results compared with the simulated response with only 
hinge danping. In the optinal control area, for the more general 
case of an asymmetrical main hub and one or two offset booms, 
the matrix Riccati equation will be solved numerically to obtain 
the control law gains. Ihe dynamics of this system will be 
simulated numerically using the optimal control law. A comparison 
in system performance both in a linear and nonlinear region will 
be made with our previous results vdiere Jhe booms were assumed to 
extend only aloig the hub principal axes. Inhere possible, for the 
. single boom offset case, a caiparison will be iiade with the previous 
results obtained by Edwards and Kaplan. 

It is planned to prepare at least one technical paper based 
on this research for open publication and/or presentation at a 
professional society meeting during the next six months. Close 
liason with NASA - Langley will be mintained. A compr^ensive 
final report will be prepared at the end of . the contractual period 
and a final oral presentation made. 

IV. POSSIBLE EXTENSION OF THE CURREJTT RESEARCH 

Referring to TABLE 1, it is hoped to conplete all the items 
related to the hinged system during the next six months. 
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Because of the success we have had with the application of linear 
control theory and the need to use the computer both for the 
solution of the matrix Ricoati equation and the simulation of the 
optinally controlled telsooping system, it is doubtful tdiat much 
progress can be made with the application of dynamic programing 
to consider the general problem of time optinal control. This is 
an item vhich could be considered in depth during a subsequent one 
year period. For all simulations involving optimal control laws it 
is assumed that such a law can be implemented with no time lag and 
no errors in measuring the state vector components - i.e. angular 
velocity conponents and boom lengths. The present work could be 
extended to include svKh an estimator problem incorporating the actual 
expected measurement errors and time lags in obtaining such measurements. 

lb date all of our results have been obtained by neglecting any 
and all external perturbations (i.e. in torque - free space). As a 
tumbling spacecraft is recovered by appropriate naneuvers of the 
appendages such effects as those due to gravity- gradient torques and 
solar pressure must be considered. It is hoped to briefly examine 
the naximum expected ixagnitudes of these perturbations during the hext 
six months. A more thoroujgh trea-ment of these two types of perturba- 
tions couild easily occupy a seconc! graduate student for a one year 
period if suitable funding were available. And finally the problem 
of incluiding flexibility effects during deployment of the telescoping 
appendages, (to the knowledge of the authors) has never been fuully treated 
in the open literature. 


This problem is made extremely difficult by the fact that the 
length of the appendage in its undefonned state is being continu£dly 
altered according to the particular control law being used. 
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TAKi: I - TOO YEAR PIAN OF STUDY 

THE DYNAMICS OF SPIN STABILIZED SPACECRAFT 
WITH MOVABLE APPENDAGES 

OONTD^TS 


A. MOTION DICING ICTLOYMENT 

Spinning spacecraft - small transverse momentum 
y/ 1. Hinged Type 

> development of equations of motion and stability analysis 
2. Telesccpic Type 

a. End mnss moving b. Uniformly distributed mass moving 

- Analytical solution for spherical Hub 

- Series solution for non- spherical Hub 
/ *3. Effect of Dampers 

B. USE OF APPENDAGES TO DETUKBLE SPACECRAFT 

1, Telescopic Type 

- derd.vation of kinetic energy 

a. Achieve zero inertial angular rate 

- Lyapunov Function - luetic Eneror 

b. Achieve spin about principal axis 

- lyapunov function - Modified kinetic energy 

4 *2. Telescoping appendages offset frcm hub principal axes 
3. Appendages + ”3" axis boom 

C. OPTIMAL CONTROL 

^ *1. Application of linear optinal 

control theory using different perfomance indices 


*2. Use of gradient technique 
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D. EFFECT or PETURBAIDNS 
*1, Gr>avity-gradient 
*2. Solar pressure 
*3. Flexibility with small amplitudle 


♦Proposed for study in second year (Part II) 




HS«23> SYSXm 


1. Derivation of Kijtetic !hergy 

2. Development of Equations of Motion 

3. IVo Dimensioned Motion Anelysis 

a. Small Angle Analysis about Equilibrium State 

(i) Linearization of Equations of Motion 

(ii) Stability Criteria 

(iii) Closed Fonn Solutions 

(iv) Nunerioal Results 

b. Lar^ Angle Analysis 

(i) Closed Fom Solutions 

(ii) Nvmerical Results 

4. Ihree Dimensioned Motion Analysis 

a. Small Angle Analysis about Equilibrium State 

(i) Linearization of Equations of Motion 

(ii) Stability Criteria 

(iii) Numerical Results 

b. Large Angle Analysis 
(i) Numerical Results 


?RECQ)ING PAGE BUNK NOT nUCQ> 


opnmi. comoL 


A. Optjjnal Control With Quadratic Perfornance Criterion 
I. Theory 

II. Applications 

1. Single Boom Offset System 

a. Equations of Mctior< 

b. Linearization 

c. Special case 

(i) Analytical Solution 

(ii) NuRi^cal Results 

- niustratlve D(aniple 

d. General Case 

(i) Nunerlcal Results 

2. TVp Boom Offset System 

a. Development of Equations of Motion 

b. Linearization 

c. Nunerioal Results 

B. Optlaal Control With Minimum Time Criterion 
I, Theory 

n. Applications 

1. Single Boom Offset System 

a. Special Case 

(i) Analytical Solution 
(11) NimeHcal Results 

b. General Case 

(i) Nmerical Results 
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